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^ ' Abstract 

h: 

. 1 We derive lower bounds for the combination of light quark masses nis + niu and 

iTid + niu- The derivation is based on first principles: the analyticity properties of two- 
point functions of local current operators and the positivity of the corresponding hadronic 
spectral functions. The bounds follow from the restriction of the sum over all possible 
hadronic states which can contribute to the hadronic spectral functions to those with the 
lowest invariant mass. These hadronic contributions are well known phenomenologically, 
either from experiment or from chiral perturbation theory calculations. The results we 
find cast serious doubts on the reliability of some of the lattice QCD determinations which 
have been recently reported. 



1 Introduction. 

There are recent determinations of the light quark masses from lattice QCD [|l|,|2| which find 
substantially lower values than those previously obtained using a variety of QCD sum rules, 
(see e.g. refs. 0-^] and earlier references therein.) The new lattice determinations are also 
in disagreement with earlier lattice results Q. More recently, an independent lattice QCD 
determination [Q] using dynamical Wilson fermions finds results which, for niu + rud are in 
agreement within errors with those of refs. [QiQli while for they agree rather well, again 
within errors, with the QCD sum rule determinations. The present situation concerning the 
values of the light quark masses is therefore rather confusing. 

The purpose of this note is to show that there exist rigorous lower bounds on how small 
the light quark masses ms+rriu and md + rriu can be. The derivation of the bounds is based on 
first principles: the analyticity properties of two-point functions of local current operators and 
the positivity of the corresponding hadronic spectral functions. Analyticity relates integrals 
of the hadronic spectral functions to the behaviour of QCD two-point functions in the deep 
euclidean region via dispersion relations. The bounds follow from the restriction of the sum 
over all possible hadronic states which can contribute to the spectral function to the state(s) 
with the lowest invariant mass. It turns out that, for the two-point functions which we 
shall consider, these hadronic contributions are well known phenomenologically; either from 
experiment or from chiral perturbation theory (xPT) calculations. On the QCD side of the 
dispersion relation, the two-point functions in question where the quark masses appear as an 
overall factor, are known in the deep euclidean region from perturbative QCD (pQCD) at the 
four loop level. The leading non-perturbative power corrections which appear in the operator 
product expansion in the physical vacuum are also known. As we shall show, the bounds 
we find cast serious doubts on the reliability of the recent lattice determinations reported in 
refs. 

We shall be concerned with two types of local operators, one is the divergence of the 
strangeness changing axial current 

^^A^^{x) = {ms + m„) :s(x)i75ii(x): , (1) 

with its corresponding two-point function 

^5(9') = i J d'xe*^ "(0|r (a^A^(x)a,^'^(0)t) |0) ; (2) 

the other one is the scalar operator S{x), defined as the isosinglet component of the mass 
term in the QCD Lagrangian 

S{x) = m[:u{x)u{x): + :d{x)d{x):], where m = -{mu + m^) , (3) 

and its corresponding two-point function 

^(^2^ = i J (i^xe^«-^(0|r (5(x)5(0)1") |0) . (4) 

We shall refer to these two-point functions respectively as the "pseudoscalar channel" and 
the "scalar channel" , and discuss them separately. 
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2 The Pseudoscalar Channel. 

The function ^'5(g^) in (Q) obeys a dispersion relation and in QCD it requires two subtractions. 
It is then appropriate to consider its second derivative (Q^ = —q^)'- 



dt 



{t + Q 



2\3 



-Im%(t), 
vr 



with 



ilm^5(t) = E(0|^M^''(0)|r)(r|5^^'^(0)^|0)(27r)3<5W(g-5:pr). 



(5) 



(6) 



The sum over T extends to ah possible on-shell hadronic states with the quantum numbers of 
the operator d^A^^ and it defines the corresponding hadronic spectral function. The behaviour 
of the ^'5(Q^) function at large euclidean values of the external momentum squared = 
—q^ ^ ^QCD known from pQCD calculations [^13|: 



87r2 Q2 



vr 



(7) 



where the dots represent higher order terms which have been calculated up to ||T3|, 
as well as non-perturbative power corrections of O (-^] |9, 14], and strange quark mass 



corrections of O and O including 0{as) terms P,p"^-[T7|. 

It is convenient to consider moments of the two-point function ^'5(9^) as follows: 



(iV + 2) 



{dQ 



N = 0,1,2, 



(8) 



and explicitly separate the contribution from the X-pole contribution to the spectral function 
in the dispersive integral in eq. (^). We then have 



(M| + Q2)iV+3 



(9) 



where 'En{Q'^) denotes the hadronic continuum integral [to = {Mk + 2m„ 



dt 



N 



to 



Im*5(t) . 



(10) 



We are now confronted with a typical moment problem |^. The positivity of the continuum 
spectral function ^lm'^^(t) constrains the moments Ti]\!{Q'^) and hence the l.h.s. of (|9|) where 
the light quark masses appear. The most general constraints among the first three moments 
for AT = 0,1,2 are: 



So(Q')>0, Si(Q2)>o, S2(Q2)>0; 
So(Q') - > 0, Si(Q2) _ S2(Q2) > ; 



(11) 
(12) 



^See e.g. ref. 
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So(Q')S2(Q2)-(Si(Q2))'>o. (13) 

The inequalities in eq. (|l^ ) are in fact trivial unless < to, which defines a region in 
where unfortunately pQCD is not applicable. The other inequalities lead however to 
interesting bounds which we next discuss. 

2a. The Bounds from T,j\f{Q'^) > 0. 
The inequality So((5^) > results in a first bound on the running masses: 

/^2n /^2x12 167r2 2/2M^ 1 1 
TUsiQ') + m.iQ^)\ >^-^U_ __„^3r^ J. (14) 



I 11 MQ^ I 
^ "I" 3 - 

The bound can of course be converted into a bound of the same combination of quark masses 
at any fixed reference scale using the renormalization group running of quark masses. The 
favourite choice of the lattice community is 4 GeV^. We shall adopt it here as well so as 
to make easier the comparison of our bounds with their results. Notice that the bound in 
eq. ( [l^ is non-trivial in the large- A'^c limit {f]^ ~ 0{Nc)) and in the chiral limit (m^ ^ ^"k)- 
The bound is of course a function of the choice of the euclidean Q-value at which the r.h.s. in 



eq. (14) is evaluated. For the bound to be meaningful, the choice of Q[GeV] has to be made 



sufficiently large so that a pQCD evaluation of the QCD factor 



-^o^"°(g') = i + ^^^ + ---, (15) 



vr 



associated to the two-point function M.\ {Q ) in eq. (^) is applicable. We consider that, 

i-Qi 




Q > 1.4 GeV is a safe choice. The evaluation of J-n'^^ at Q = 1.4 GeV with or without 



inclusion of the 0{al) terms differs by less than 10%. The lower bound which follows from 
eq. (0) for ruu+ms at a renormalization scale /U^ = 4GeV'^ results in the solid curves shown in 
Fig. 1 below. These are the lower bounds obtained by letting vary [19[] between 290 MeV 
(the upper curve) and 380 Mev (the lower curve) and using fx = 113MeV, Mk = 493.67MeV. 
We take the number of active flavours Uf = 3 in all the numerical analyses of this work. Notice 
that the hadronic continuum integral in ( p!o[ ) is always larger than the contribution from three 
light flavours; and indeed, we have checked that using Uf = 4 in the QCD expressions leads to 
even higher bounds. Values of the quark masses below the solid curves in Fig. 1 are forbidden 
by the bounds. For comparison, the horizontal lines in the figure correspond to the central 
values obtained by the authors of ref. Q: their "quenched" result is the horizontal upper line; 
their "unquenched" result the horizontal lower line. 

It should be emphasized that in Fig. 1 we are comparing what is meant to be a "calcula- 
tion" of the quark masses -the horizontal lines which are the lattice results of ref. 0- with 
a bound which in fact can only be saturated in the limit where So(Q^) = 0. Physically, this 
limit corresponds to the extreme case where the hadronic spectral function from the contin- 
uum of states is totally neglected! What the plots show is that, even in that extreme limit, 
and for values of Q in the range 1.4 GeV ^ Q ^ 1.8 GeV the new lattice results are already 
in difficulty with the bounds. 

■^Lower bounds on light quark masses based on inequalities of this type were first discussed in ref. Q. Our 
discussion here is an update of this work. 
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Fig. 1 Lower bound in MeV for [mg + Tnu]{4: GeV ) versus Q(GeV) /rom eg. fi^j. 



The inequalities > and T,2{Q'^) > result in the improved bounds 



2 ^ levr^ 2/2.M^ 



1 



1 + ^ 



1 I 17 Q:4Q^) I 



(16) 



and 



2 167r2 2/2 M|. 6 



> 



1 + 4^ 



1 I 20 MQ^ I 
^ 3 vr 



(17) 



The improvement is due to the fact that higher and higher moments increase the weight of 
the contribution from the low energy part of the hadronic spectrum, in our case the ii'-pole; 
e.g., at Q = 2 GeV and for A^^^ = 330 MeV eq.(|g) imphes (m^ + m„)[4 GeV^] > 77 MeV 



and eq.(17) (m^ + m„)[4 GeV^] > 106 MeV. Notice however that as the bounds are improved 
by going to higher moments, the corresponding QCD functions 



1 I 17 «s(Q^) ^ 



^2''^°(Q2) 



^ I 20 a.(Q2) ^ ^ 



vr 



(18) 



3 vr z v-v . 3 

have larger coefficients in their pQCD expansion in powers of ag. Therefore, for the improved 
bounds to be meaningful, the appropriate choices of Q have to be made at sufficiently larger 
values so that a pQCD evaluation of the J-'^*^°-functions is applicable. The same phenomenon 
occurs when one considers the quadratic bound in eq.(|l3|) which we next discuss. However, as 
we shall see, the net improvement at the appropriate higher Q-values is still quite considerable. 
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2a. The Bounds from the Quadratic Inequality. 
Inserting eq. (P) in the quadratic inequality in eq. (|^) and solving for the quark masses 
results in the bound 

16^2 



rUsiQ ) + mu{Q 



> 



Nr. 



9 

A/ 2 



■QCD 
2 



(19) 



the two-point functions A^5^^ in eq.(S) with = 0, 1, 2: 



where -Fo'^^°(Q2), jrQCD(Q2-)^ jr2«^D(Q2^ are the QCD factors (|l|) and (H) associated to 

5 

QCD^2 _ ^ , , ^20) 



q T^QCD ^QCD 
"J-^ -^2 



o vr 



The resulting lower bounds for the quark masses are the solid lines shown in Fig. 2 below. 
Like in Fig. 1, these are the lower bounds obtained by letting vary [l^ between 290 MeV 
(the upper curve) and 380 Mev (the lower curve). Values of the quark masses below these 
solid curves are forbidden by the bounds. The horizontal lines in this figure are also the same 
lattice results as in Fig. 1. 
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Q[GeV] 



Fig. 2 Lower bound in MeV for [mg + (4 GeV ) from the quadratic inequality in eq. ([Z 



The quadratic bound in eq. (jl|) is saturated when So(Q2)S2(Q2) _ {T.i{Q'^)) = 0, which 



happens for a 5-like spectral function representation of the hadronic continuum of states at 
an arbitrary position and with an arbitrary weight. This is certainly less restrictive than the 
extreme limit with the full hadronic continuum neglected, and it is therefore not surprising 
that the quadratic bound happens to be better than the ones from S7v(Q^) for = 0, 1, and 
2. Notice that the quadratic bound in Fig. 2 is plotted at higher Q- values than the bound in 
Fig. 1. The evaluation of the QCD factor in (pO|) with or without inclusion of the 0(af ) terms 
differs by less than 10% for > 4 QeV^ and we consider therefore that for the evaluation 
of the quadratic bound Q > 2 GeV is already a safe choice. We find that even the quenched 
lattice results of refs. are in difficulty with these bounds. 

Similar bounds can be obtained for m„ + m^^ if instead of the two-point function in (|2|) 
one considers the two-point function associated to the divergence of the axial current 

^^A^'{x) = {rud + m„) ■.d{x)i-i^u{x): . (21) 

The method to derive the bounds is exactly the same as the one discussed above and therefore 
we only show, in Fig. 3 below, the results for the corresponding lower bounds which we obtain 
from the quadratic inequality. 




We find again that the lattice QCD results of refs. [Q,^,^] for m^ + m^ are in serious difficulties 
with these bounds. Notice that in this case the mass corrections in the QCD two-point 
function are negligible at the plotted Q-values. 
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3 The Scalar Channel. 



The function 'if{q'^) in eq. (Q) obeys a dispersion relation which in QCD requires two subtrac- 
tions. We shall thus work with its second derivative 



dt- 



1 



-lm^(t) . 



(22) 



The QCD behaviour of the "^"(Q^) function at large euclidean values of the external mo- 
mentum squared = —q^ ^ ^Qcr> entirely similar to the one of '^'^^{Q'^) which we have 
already discussed, but for very small m„ and quark mass correction terms which we now 
neglect i.e., 



2 X 



87r2 Q2 



1 + 



llas(Q' 



vr 



+ ... 



(23) 



The difference between the "scalar" and "pseudoscalar" channels is that the lowest state 
which contributes to the hadronic spectral function in the "scalar" channel is not a pole but 
the vr — vr continuum with J = and / = 0. This contribution provides a lower bound to the 
full spectral function: 



^3|F(i)|Mi-4m^), 



where F{t) denotes the scalar-isoscalar pion form factor 

{'K\p)TT\p')\S{Qm = 6abF{t) , With 



t = {jp+p' 



l\2 



(24) 



(25) 



In particular, the value of this form factor at the origin is the so called pion sigma term. This 
form factor has been well studied in xPT at the one and two loop level (see e.g. refs 21 
and references therein), and the predicted low energy shape agrees well with the available 
experimental information. With 



Fit) = F{0) 



1 + ^{r^^t + 0{t'] 



it is found M that 



F{0) 



X 



(0.99 ± 0.02) + 0{ml) and (r^)J = (0.59 ± 0.02) fm 



(26) 



(27) 



Using standard methods ^ we construct bounds on m when something is known about the 
scalar-isoscalar pion form factor F(t). For this, it is convenient to map the complex g^-plane 
to the complex unit disc: 



'l-z 



with 



Ami 



(28) 



^The pion mass used in ref. is the 7r^-mass. Since the electromagnetic interactions are neglected, we 
find it more appropriate to use the vr^-mass instead; which in any case results in a lower contribution to the 
average quark mass. In the generalized version of xPT (see e.g. ref. and references therein,) the value of 
F{0) could be sensibly larger, which would result in even larger bounds. A value for the curvature of the scalar 
form factor is also quoted in ref. ; it involves however some extra phenomenological assumptions which go 
beyond xPT and this is why we have not included this extra information here. 

■^see ref. |^ and references therein, and also refs. [p^p6[ 
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The dispersion relation in (p2|) and the spectral function inequality in ( p4| ) can then be 
combined as follows: 



W{Q')>lj^dy[y + -^^ y-"Hy-lfl'\F{Amly)\^^ J ^ ^|,^(e-)F(^(e-))|^ 

(29) 

where (piz) is an analytic function on the unit disc which contains the phase-space factors as 
well as the dispersive weight factor and the jacobian of the conformal mapping; and 



levr 



167r2 Q2 



vr 



(30) 



With F(0) and F'{0) as input, the analiticity of (f){z)F{z) on the unit disc enables us to 
write the constraint |25 



>V(Q2) > 



1 + 



.F\0) + 



4ml) 

which results in a scalar channel bound for iv?: 



(31) 



^2r^2.^ ^ rnl [AmlY 3 x (0.99 ± 0.02)^ I + {-^ + H-Q'] - '-^{rYsY 



1 I 11 I 

3 TT 



Notice that this bound is somewhat similar to the first bound obtained from the divergence 
of the axial current in eq. ( [I^ ) . There are however some interesting differences which we wish 
to point out. From the point of view of the large- A^c expansion, the r.h.s. of eq. ( |3^ ) is 1/Nc- 
suppressed. Also, from the point of view of the chiral expansion the r.h.s. of eq. ( |3^ has an 
extra factor m"^ as compared to the r.h.s. of eq. ([l^. One might expect therefore that the 
lower bound in (32) should be much worse than the one in (IH 



Yet the bound for m in eq.(p2[) 
is surprisingly competitive. The reason for this is twofold. On the one hand, the absence of 
a prominent narrow resonance in the J = 0, / = channel, is compensated by a rather large 
contribution from the chiral loops, a phenomenological feature which appears to be common 
in all the explored physical processes where the J = 0, / = channel contributes. On the 
other hand the fact that, numerically, and are not so different. The lower bounds 
for the sum of running masses (m„ + m^) at the renormalization scale [/i^ = 4 GeV^] which 
follow from this scalar channel bound are the dashed curves shown in Fig 4, corresponding 
to the choice = 290 MeV (the upper dashed curve) and A^ = 380 MeV (the lower 

iVi J Ivl O 

dashed curve). The plots show that for values of Q in the range 1.4 GeV Q ^ 1-6 GeV the 
lattice QCD results in |jl],0,0 are dangerously close to lower bounds obtained with very little 
phenomenological input. Notice also that the relevant two-point function in this channel is 
practically the same as Ai^^\Q^) in eq. (|8|) for which Q > 1.4 GeV is already a safe choice. 

One can still improve these bounds by taking into account the extra information provided 
by the fact that the phase of the scalar form factor F{t) in the elastic region 4m2 < t < 16m2 
is precisely the / = 0, J = tt — tt scattering phase-shift and there is information on this 
phase-shift both from xPT and from experiment. 
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1.4 1.6 1.8 2 2.2 2.4 

Q[GeV] 

Fig. 4 Lower bounds in MeV /or [m^ + mu](4 GeV^) /rom the Scalar Channel. 



The technology to incorporate this information is discussed in refs. |24-Pq]. We have restricted 
the phase-shift input to a conservative region: 4m^ < i < 500 MeV^, where a resummation of 
the chiral logarithms calculated at the two-loop level in xPT [^,^ is certainly expected to 
be reliable and is in good agreement with experiment. This results in improved bounds which 
correspond to the solid curves shown in Fig. 4. The two curves reflect the variation within 
the quoted errors of the input parameters. Again, we find that the lattice determinations of 
refs. Il|,|2l0] do not satisfy lower bounds which only incorporate very little hadronic input. 



4 Conclusions. 

We conclude that the lattice results of refs. are in serious difficulties with rigorous 

lower bounds which can be derived from general properties of QCD and with a minimum of 
well established phenomenological input; and this in two different hadronic channels. The 
lower bounds we have derived are perfectly compatible with the sum rules results of refs. |^ 
^] and earlier references therein. They are also compatible with a recent semiempirical 
determination of made by the ALEPH collaboration |^] based on the rate of the hadronic 
r-decays into strange particles observed at LEP. 

^There exists a recent preprint |Q with various criticisms on the determination of quark masses with QCD 
sum rules. A reply to these criticisms will be made elsewhere |28[. 
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